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Abstract. The nuclei of active galaxies harbor massive young stars, an accreting central black hole, or both. In order to deter- 
mine the physical conditions that pertain to molecular gas close to the sources of radiation, numerical models are constructed. 
These models iteratively determine the thermal and chemical balance of molecular gas that is exposed to X-rays (1-100 keV) 
and far-ultraviolet radiation (6-13.6 eV), as a function of depth. We present a grid of XDR and PDR models that span ranges 
in density (10 2 - 10 6 ' 5 cm" 3 ), irradiation (10 05 - 10 5 G and Fx = 1.6 x 10" 2 - 160 erg cm" 2 s" 1 ) and column density 
(3 x 10 21 — 1 x 10 25 cm -2 ). Predictions are made for the most important atomic fine- structure lines, e.g., [CII], [OI], [CI], 
[Sill], and for molecular species like HCO+, HCN, HNC, CS and SiO up to J = 4, CO and 13 CO up to J = 16, and col- 
umn densities for CN, CH, CH + , HCO, HOC + , NO and N2H + . We find that surface temperatures are higher (lower) in PDRs 
compared to XDRs for densities > 10 4 (< 10 4 ) cm" 3 . For the atomic lines, we find that, largely due to the different XDR 
ionization balance, the fine- structure line ratios of [Sill] 35 ^im/[CII] 158 /xm, [OI] 63 ^tm/[CII] 158 fim, [Fell] 26 /im/[CII] 
158 /im and [CI] 369 /im/[CI] 609 /im are larger in XDRs than in PDRs, for a given density, column and irradiation strength. 
Similarly, for the molecular lines, we find that the line ratios HCN/HCO + and HNC/HCN, as well as the column density ratio 
CN/HCN, discriminate between PDRs and XDRs. In particular, the HCN/HCO+ 1-0 ratio is < 1 (> 1) for XDRs (PDRs) if the 
density exceeds 10 5 cm" 3 and if the column density is larger than 10 23 cm" 2 . For columns less than 10 22 ' 5 cm" 2 the XDR 
HCN/HCO + 1-0 ratio becomes larger than one, although the individual HCN 1-0 and HCO + 1-0 line intensities are weaker. 
For modest densities, n = 10 4 — 10 5 cm" 3 , and strong radiation fields (> 100 erg s" 1 cm" 2 ), HCN/HCO + ratios can become 
larger in XDRs than PDRs as well. Also, the HCN/CO 1-0 ratio is typically smaller in XDRs, and the HCN emission in XDRs is 
boosted with respect to CO only for high (column) density gas, with columns in excess of 10 23 cm" 2 and densities larger than 
10 4 cm" 3 . Furthermore, CO is typically warmer in XDRs than in PDRs, for the same total energy input. This leads to higher 
CO J=N+l-N/CO 1-0, N>1, line ratios in XDRs. In particular, lines with TV > 10, like CO(16-15) and CO(10-9) observable 
with HIFI/Herschel, discriminate very well between XDRs and PDRs. This is crucial since the XDR/AGN contribution will 
typically be of a much smaller (possibly beam diluted) angular scale and a 10-25% PDR contribution can already suppress 
XDR distinguishing features involving HCN/HCO+ and HNC/HCN. For possible future observations, column density ratios 
indicate that CH, CH+, NO, HOC+ and HCO are good PDR/XDR discriminators. 



1. Introduction 

The radiation that emanates from galaxy nuclei, such as those 
of NGC 253 and NGC 1068 or even more extreme (ultra- 
luminous infrared galaxies, is believed to originate from re- 
gions with active star formation, an accretion disk around 
a central super -massive black hole or bo th (e.g. Silk |2005; 
Maloney 1999; Sanders & Mirabel 1996). The unambiguous 
identification of the central energy source, or the relative con- 
tributions from stars and an active galaxy nucleus (AGN), re- 
mains a major challenge in the study of active galaxy centers. 
The general aim of this paper is to determine how the proper- 
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ties of the irradiated interstellar medium (ISM) may further our 
understanding. 

Unlike emission at optical wavelengths, atomic, molecular 
and dust emission in the far-infrared and (sub-)mm range allow 
one to probe deeply into the large column densities of gas and 
dust that occupy the centers of active galaxies. Observational 
studies of the ISM in galaxy centers h ave be en present ed by 
vario us auth ors (c.f. Aalto 2005; BaanEoOS Ott et al.|2005; 
Israel l2005t Spoon et al. 1200 it l2003l bo05t K lockner et al. 
2003| Isr ael and Baas l2002t Garret t et al. 1200 it Hiittemeister 
and Aalto 2001; Curran et al. 2000). Theoretical models show 
that the spectral energy distribution of the radiation represent- 
ing star formation (peaking in the ultraviolet) and AGN (peak- 
ing in the X-ray regime) activity respectively influences the 
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thermal and chemical balance of the ambient ISM in funda- 
mentally different ways (Meijerink & Spaans 2005; Maloney, 
Hollenbach & Tielens 1996; Lepp & Dalgarno 1996). The spe- 
cific aim of this paper is thus to study the extent to which emis- 
sion from commonly observed molecular and atomic line tran- 
sitions may be used as a diagnostic tool in the study of external 
galaxy centers to determine the ambient conditions in general, 
and the type of irradiation in particular. 

To this effect, we have extended the chemical calculations 
described by Meijerink & Spaans (2005; hereafter Paper I) for 
ultraviolet and X-ray irradiated gas to a much larger parameter 
space of ambient conditions and we have performed detailed 
radiative transfer calculations to compute the line intensities 
of many atomic and molecular transitions. We refer the inter- 
ested reader to Paper I for a detailed description of the com- 
bined photon-dominated region (PDR)/X-ray dominated region 
(XDR) code that we have used to compute the impact of ultra- 
violet (PDR) and X-ray (XDR) photons on nearby ISM. The 
results described here will be applied to observations of ex- 
ternal galaxy c enters in a subsequen t paper (Meijerink et al., 
in preparation). Staube r~t alJ J2005h also developed PDR and 
XDR codes, for applications to Young Stellar Objects. The 
main difference is that this code has a time dependent chem- 
istry and that it includes ice evaporation at t = 0. This affects 
the chemistry compared with traditional steady-state models of 
pure gas-phase PDRs (e.g., CN/HCN). This is not relevant for 
our cloud models. However, freeze-out and evaporation start to 
become important for clouds with density na > 10 5 cm -3 . 

2. 'Standard' clouds 

The current spatial resolution of sub-millimeter telescopes, 
such as the James Clerk Maxwell Telescope (JCMT), the 
Institute de Radio Astronomie Millimetrique (IRAM) telescope 
and even the Combined Array for Research in Submillimeter 
Astronomy (CARMA) or the Submillimeter Array (SMA), 
is insufficient to resolve individual clouds in extragalactic 
sources. By using these telescopes, each resolution element 
thus measures the combined emission from a large ensemble of 
molecular clouds. As a consequence, it is frequently impossible 
to use a single model cloud solution to describe the observed 
molecular lines. Instead, more complicated solutions involving 
two or more model clouds, with differing densities and inci- 
dent radiation fields, are needed. This contrasts with the study 
of Galactic objects, where usually a single model cloud solu- 
tion is sufficient to fit the measurements of single resolution 
elements. 

In this paper, we calculate a grid of 'standard' clouds sam- 
pling the different physical conditions believed to be relevant 
for the centers of active galaxy nuclei. In order to sample both 
the hierarchical size and (column) density properties of the 
ISM, we have chosen to construct models for a number of fixed 
sizes as well as densities. Note that the column densities are 
not the same for each model, since we use fixed cloud sizes. In 
these clouds, we investigate the detailed column density depen- 
dence for the line ratios HCN/CO, HNC/HCN, HCO+/HCN, 
SiO/CO and CS/HCN, which include line ratios observed in 
several galaxies. From our computational grid, one or more, 



properly weighted, models can be chosen to reproduce ob- 
served atomic and molecular lines. 

We distinguish three different 'standard' clouds, each with 
their own characteristic combinations of size and volume 
density range, hence also column density (cf. Table [!}, for 
which we calculate a set of models for different incident 
radiation fields, and where a distinction between irradiation 
by far-ultraviolet (FUV) and X-ray photons is made. The 
X-ray radiation field is a power-law distribution F(E) = 
F(0)(E/lkeV)~ a integrated between 1 and 100 keV, where 
a = —0.9. A low energy cut-off of 1 keV is chosen, since 
most of the emitted radiation below 1 keV is absorbed in the 
hot, highly ionized medium close to the black hole. This power 
law spectrum is generally believed to be representative for ac- 
creting black holes. Please note that this differs from Paper I, 
where we used a thermal distribution at 10 6 K instead. The ul- 
traviolet radiation field (6-13.6 eV) is expressed in multiples 
of the Habing flux, where Go — 1 corresponds to 1.6 x 10 -3 
erg cm -2 s _1 , which is the local Galactic interstellar radiation 
field. We express the X-ray flux in erg s _1 cm -2 . An accreting 
black hole with an X-ray luminosity of 10 44 erg s _1 , produces 
an X-ray flux of 100 erg -1 cm -2 at a distance of 100 pc, when 
there is not extinction. We use a line width 5 V — 2.7 km s _1 . 
Cloud type A represents compact, high-density environments 
such as molecular cloud cores, and clouds very close to active 
nuclei; cloud type B corresponds to more traditional molecular 
cloud environments, and cloud type C is representative of the 
more diffuse extended (molecular) medium in which clouds of 
type B are usually embedded. 

Late-type galaxies frequently have radial metallic- 
ity gradients, with the highest metallicity in th e center 
dVila-Costas & Edmunds! 1 1992HZaritskv et all 1 1994 . For this 
reason, we have adjusted the metallicity used in Table 2 of 
Paper I. The published metallicity gradients and suggestions 
of a gradient flattening in the very center have led us to adopt 
a twice Solar metallicity as a reasonable value. Since [C]/[0] 
abundance ratios decrease at higher metallicities, we have 
taken the carbon abundance e qual to the oxygen abundance 
(see, for instance, Garnett et al. 1999; Kobulnicky & Skillman 
1998). Note that the [Cl/[ qi ratio affects the abundances o f 2 
and H2O. See for e xample | $paans ffi: van Dishoeckl fcoOl) and 
especially Fig. 2 in Bergin et al. (2000). 

From the models, we have calculated the intensities of the 
molecular rotational lines of HCN, HNC, HCO+, CS and SiO 
(upto J = 4), [CI], [CII], [OI], [Sill] and other fine-structure 
lines. For CO and 13 CO we calculated the intensities of the ro- 
tational lines up to J = 16, in order to make predictions for fu- 
ture observing facilities such as the ESO Herschel/HIFI space 
mission. We use the Leid en Atomic and Mole cular Database 
(LAMDA) as described in Schoi er et ail J2005I) to retrieve the 
collisional data needed for the calculations. Where no colli- 
sional data are available for commonly observed molecules 
such as CN, CH, CH+, HCO, HOC+, NO and N 2 H+, we only 
give column densities. 
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Type 


r(pc) 


7i(cm 
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iV(cm- 2 
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Go 




F x [erg s 


cm" 2 ] 


A 


1 


10 4 - 


1Q 6.5 


3 x 10 22 


- 1 x 10 25 


10 2 - 


10 5 


1.6-160 




B 


10 


10 3 - 


10 4 


3 x 10 22 


- 3 x 10 23 


10 1 - 


10 4 


1.6 x 10" 


x -16 


C 


10 


10 2 - 


10 3 


3 x 10 21 


- 3 x 10 22 


1Q 0.5 


- 10 3 


1.6 x 10" 


2 -1.6 



3. Surface temperatures 

To illustrate the coupling differences of FUV and X-ray pho- 
tons to the gas, we first discuss the surface temperatures of the 
low, mid, and high density models respectively over the param- 
eter space given for the PDR and XDR models in Tabled We 
calculated a larger range of radiation fields for the PDR mod- 
els than for the XDR models, for reasons related to the heat- 
ing efficiency which are discussed below. The resulting surface 
temperatures as a function of gas density and incident radiation 
intensity for both PDR and XDR models are shown in Fig.^ 

The most important heating mechanism at the edge of a 
PDR is photo-electric heating. FUV photons are absorbed by 
dust grains, which release electrons that lose their surplus ki- 
netic energy to the gas by Coulomb interactions. The efficiency 
of this process increases, when the grains are more negatively 
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Fig. 1. Surface temperatures for PDR (top) and XDR (bottom) 
models. 



charged, which is determined by a complex interplay of the im- 
pinging radiation field Go, electron density n e and gas temper- 
ature T. The absolute efficiency is very low, since about 0.5-3 
percent of the photon energy goes into gas heating. This sharply 
contrasts with direct X-ray heating, important at the edge of 
XDRs. Direct ionization of an atom yields a kinetic electron 
with an energy higher than 1 keV. This electron heats, ionizes 
and excites the gas. Depending on the H, H2, He and elec- 
tron abundances, the heating efficiency can be up to 70 percent, 
much higher than for photo-electric heating. However, there is 
an opposing effect, namely the much lower absorption cross 
section for X-rays. Since the cross sections scale roughly as 
E~ 3 , there are many fewer X-ray photons absorbed than FUV 
photons. 

We find that at high densities (n > 10 4 cm -3 ) PDR models 
produce higher surface temperatures than the XDR models. At 
low densities, however, we find the opposite, especially in the 
case of high radiation fields. This is explained by the drop in 
the efficiency of photo-electric heating due to grain charging 
at densities n < 10 5 cm -3 when the same impinging radia- 
tion field is considere d. This is nicely illustrated in Fig. 6 of 
iKaufman et alJ lll999l) . where the ratios of the intensity of the 
[CII] 158 /im and [OI] 63 /im lines to the total far-infrared in- 
tensity emitted from the surface of the clouds are plotted as a 
function of density and radiation field. This ratio is a measure 
of the heating efficiency, since [CII] 158 /im and [OI] 63 /im 
are the most important coolants in PDRs. 

In the regime discussed here, the surface temperatures of 
the XDR models are quite well correlated with Hx I n, where 
Hx is the energy deposition per particle and n the total hy- 
drogen density. This means that the highest surface temper- 
atures are found at the lowest densities and highest imping- 
ing radiation fields. Consequently, the contours of equal sur- 
face temperature are almost straight lines in the XDR plot 
of Fig. [l] In the PDR models their beh avior is more com- 
plicated, as already discussed bv lKaufman et alJ \ 19991) . since 
grains are involved in heating the gas. At the edge of the cloud, 
the cooling is dominated by [CII] 158 and [OI] 63 /im, 
which have critical densities of n cr (CU) rj 3 x 10 3 cm~ 3 and 
n cr (OI) « 5 x 10 5 cm -3 for collisions with atomic and molec- 
ular hydrogen. The cooling rate in this regime is more or less 
proportional to n 2 . The heating rate is at least proportional to n, 
because the grain density is proportional to n. It can be larger, 
because grains become less positively charged at increasing 
electron densities making it easier for electrons to escape the 
grains, so that the heating efficiency increases. Nevertheless, 
the density dependence of heating remains less steep than n 2 , 
i.e., less than that of cooling, which causes a drop in the tem- 
perature at n < 10 4 cm~ 3 and for fixed Go with increasing 
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density. Between n m 1—3 x 10 3 cm" 3 where the [CII] 158 /im 
line thermalizes, the drop in temperature stagnates. For densi- 
ties between n = 10 4 — 10 5 cm -3 and Go < 10 4 , we also find 
that the surface temperature drops with increasing n. The [OI] 
63/im line thermalizes at higher densities n — 10 5 — 10 6 cm -3 
and, therefore, cooling will be proportional to n. Heating in- 
creases faster with density, which results in a higher surface 
temperature. For Go > 10 4 , we find that the surface tempera- 
ture rises up to n ~ 5 x 10 5 cm -3 due to the increase in the 
heating efficiency with density at fixed Go- Above this density 
and at these high temperatures, coolants with high critical den- 
sities and excitation energies such as the [OI] 6300A line be- 
come important, causing the surface temperature to drop again. 

Wh en we compare the surf ace temperatures with those de- 
rived by Kauf man et alJ Jl999l) . we find that our model surface 
temperatures are higher. A possi ble exp lanation for this has 
been discussed recently by Roll ig et alJ (|2006), who present 
scaling relations for heating and cooling as a function of metal- 
licity Z. They state that the photo-electric heating rate is oc Z 
for n < 10 3 cm -3 increasing to Z 2 when n > 10 6 cm" 3 . The 
cooling rate is always proportional to Z, and, therefore, higher 
metallicities result in higher surface temperatures. Note how- 
ever that the temperature differences found are very likely not 
only bec ause of a change in metallicity. In the PDR compar- 
ison test jRollig et all 12006). we found significant scatter be- 
tween different PDR codes in the thermal balance. Therefore, 
one should not take the absolute values of line intensities too 
literally in the interpretation of data. 

4. Fine-structure lines 

In principle, we can use combinations of fine-structure line 
intensities to constrain densities and incident radiation fields. 
Here we discuss how such line ratios depend on gas density, 
ambient radiation field strength, and cloud column density, and 
compare the results for the PDR and XDR mo dels. In previ- 
ous st udies only results for PDRs were shown. IWolfire et alJ 
Jl990h were the first to show r CII]/[ SiII1 contour plots for 
PDRs. Recently, Kauf man et alJ ll2006h showed updated plots 
for [Sill] and [Fell] for PDR gas. The fine-structure lines are 
selfconsistently calculated in the PDR and XDR codes, and we 
take line-trapping into account for ALL lines. 

4. 1. [Sill] 35 fj,m/[CII] 158 fim intensity ratio 

On the website, we show the [Sill] 34.8 /im / [CII] 158 /im fine- 
structure line ratio for both the PDR and the XDR models. The 
[Sill] 34.8 /im line has an energy of E/k = 414 K and a criti- 
cal density of n cr = 3.4 x 10 5 cm" 3 , while this is E/k = 92 K 
and n cr = 2.8 x 10 3 cm" 3 for the [CII] 158 /im. Very high ra- 
diation fields produce gas temperatures in the PDR models that 
are sufficiently high to excite the [Sill] and [CII] upper states 
and the ratio depends mostly on density. With lower radiation 
fields, the surface temperature drops. The upper level energy 
of [Sill] 34.8 /im is reached first and the [SiII]/[CII] line ratio 
drops. As the density increases at a given FUV radiation field, 
the ratios limit to roughly constant values that are set by the 
corresponding surface temperatures. 



In the XDR models, the ratio is not only determined by the 
temperature and density, but also by the fractional abundances 
and column densities. In PDRs clearly defined layers occur, 
in which carbon and silicon are both almost fully ionized, but 
throughout XDRs neutral and ionized species co-exist. Despite 
the fact that surface temperatures in the XDR models are lower, 
we find much higher [SiII]/[CII] ratios, since silicon is much 
easier to ionize than carbon. The dominant source for ioniza- 
tion is not the direct absorption of an X-ray photon (primary 
ionization), but the produced kinetic electron. This electron can 
ionize a species either directly by collisions (secondary ioniza- 
tions) or indirectly by first exciting H and H2 and producing 
Lyman a and Lyman- Werner photons, which then may ionize 
species in turn. The cross section of Si for secondary ioniza- 
tion is about twice that of C. This does not, however, fully ex- 
plain the order of magnitude difference between the calculated 
PDR and XDR ratios. This difference also reflects the fact that 
ionization of C can only be done by Lyman- Werner photons, 
whereas both Lyman a and Lyman- Werner photons are capa- 
ble of ionizing Si. It is thus harder to ionize C than Si in regions 
where the gas is mostly atomic, which results in a large increase 
of the ratio for all densities at a given irradiation strength. 

This also explains the fairly constant ratio below Fx ~ 
10 erg s" 1 cm" 2 for n ~ 2 x 10 3 — 3 x 10 6 cm" 3 , which re- 
sults from an interplay between various effects. With ambient 
radiation fields constant, we find that the surface temperatures 
drop with increasing density in the XDR. At lower tempera- 
tures and higher densities, H2 is more easily formed. Both the 
temperature drop itself and the enhanced H2 (leading to more 
carbon ionizations) thus suppresses the [SiII]/[CII] ratio, but 
this is compensated by the relatively high critical density of the 
[SiH] 34.8 /xm line. 

For the highest, > 16 erg s" 1 cm -2 , radiation fields, where 
carbon is highly ionized, we find the same trend as seen for 
the PDR, i.e., the ratio is mostly dependent on density for the 
highest radiation fields, and then there is a decrease in the ratio 
when the temperature drops below the upper-state energy of 
[Sill]. 

At the lowest densities and highest X-ray radiation fields, 
we find that the effect of column density become important as 
well. Since we fixed cloud sizes, in each cloud type the lowest 
density models imply also the lowest column densities. In the 
high irradiation models, carbon is almost fully ionized at the 
XDR edge. The fractional abundance of ionized carbon drops 
toward the H/H2 transition and then increases again for the rea- 
sons discussed earlier. In the lowest (column) density models, 
we only produce the highly ionized part, suppressing the ra- 
tio even more. This is best seen in the diagram for the high 
density XDR models, at densities n = 10 4 cm" 3 . The ratio 
increases from Fx = L6 to 30 erg s" 1 cm" 2 and then drops 
again. Hence, if we increase the column densities of these mod- 
els by a factor of ~ 30, the ratios only depend on radiation field 
strength and not on density. 
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Fig. 2. [OI] 63 /xm / [CII] 158 /xm ratio for PDR (left) and XDR (right) models. 



4.2. /O/; 63 fim/[CII] 158 urn intensity ratio 

In Fig. El we show the [OI] 63 /xm / [CII] 158 jtxm ratio. The 
critical density 10 cm and the upper-state en- 

ergy is E/k = 228 K for [OI] 63 /xm. In the PDR models and 
at very high incident radiation fields, this ratio depends mostly 
on density as in the case for the [SiII]/[CII] ratio, since once 
again temperatures are sufficiently high to excite both upper- 
state levels. When the surface temperature drops for lower ra- 
diation fields, the upper-state energy of [OI] is reached first. 
The [OI]/[CII] ratio then drops, with a flat density dependence 
for n < 10 4 cm~ 3 and Go < 10 3 . The decrease in this line ra- 
tio is, however, less pronounced than that in the [SiII]/[CII] ra- 
tio, because the upper-state energy of [OI] is lower than that of 
[Sill]. When G < 10 3 and rx> 10 5 5 cm~ 3 , the ratio is roughly 
constant at the same density for all values of Go- Here, the 
more difficult excitation of [OI] at lower temperatures is coun- 
teracted by the C + layer becoming thinner at lower radiation 
fields. 



In the XDR, the [OI]/[CII] ratio is again more complex 
due to its dependence on density, radiation field, ionized car- 
bon fraction and total column density. The ratios are overall 
much higher than in the PDR models, since carbon does not 
become fully ionized. At the largest Hx /n, temperatures are 
high enough to create a dependence on density only. At lower 
radiation fields, the [OI]/[CII] ratio drops as temperatures ap- 
proach the upper-state energy of the [OI] 63 fim line. When 
we decrease the radiation field even more, we find that the C + 
fraction is rapidly reduced at high (> 10 5 cm~ 3 ) densities and 
the ratio increases again. At high (> 16 erg s _1 cm~ 2 ) radi- 
ation field strength and low densities, a large Hx / n is main- 
tained throughout the whole cloud, since the relevant type A 
cloud size is only one parsec. For that reason, we find here 
the same effect as already seen in the [SiII]/[CII] ratio, since 
column densities increase toward higher densities. At the high- 
est radiation fields and lowest densities, the ratio is suppressed 
since carbon remains partially ionized over the full extent of 
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Low density PDRs: [CI] 569fim / [Cl] 609/itm 



Low density XDRs: [CI] 369/im / [CI] 609/im 
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Fig. 3. [CI] 369 fim I [CI] 609 fim ratio for PDR (left) and XDR (right) models. 



the clouds considered here. We do not find the highest ratio at 
the highest radiation field in the lower density models. 

4.3. [Fell] 26 fim/[CII] 158 fim intensity ratio 

On the website, we show the [Fell] 26 fim I [CII] 158 fim in- 
tensity ratio. [Fell] 26 fim is very difficult to excite due to its 
high critical density n cr = 2.2 x 10 6 cm" 3 and upper-state en- 
ergy E/k — 554 K. Thus, the change in ratio with increasing 
density at high incident radiation fields is much larger than that 
seen for the [Sill] 34.8 fim I [CII] 158 ^m ratio. In the PDR 
models and at high incident radiation fields, the ratio mostly 
depends on the density. At lower radiation field strengths, ap- 
proaching the upper-state energy of [Fell], the ratio drops. 

The same trends are seen at high radiation fields in the XDR 
models, but again, we find much higher ratios than in the PDR 
models. It is possible to ionize iron with Lyman a photons, but 
not carbon. At moderate radiation fields, we again find ratios to 
be more or less independent of density, for the reason that we 



have already discussed in the [SiII]/[CII] case. At the lowest 
densities and highest radiation fields in each cloud type, we 
find that the ratios are similarly suppressed as was the case for 
[Sill] 35 /Ltm/[CII] 158 fim and [OI] 63 /zm/[CII] 158 fim. 

4.4. [Cl] 369 fim/[CI] 609 fim intensity ratio 

In Fig.|3j we show the [CI] 369 fim I [CI] 609 fim intensity ra- 
tio. The critical densities of these lines are n cr = 3 x 10 2 cm -2 
for [CI] 609/rni and n cr = 2 x 10 3 cm" 2 for [CI] 369 fim, 
typically lower than the densities we are considering here. The 
upper-state energies are E/k = 24 and 63 K for [CI] 609/im 
and [CI] 369 fim, respectively. In the PDR models, the [CI] 
lines originate from the C+/C/CO transition layer. The temper- 
atures in this layer slightly rise with increasing incident radia- 
tion field strengths and range between T = 20 — 100 K (com- 
parable to the upper-state energies). This explains the small in- 
crease in the ratio for larger Go at the same density. The ra- 
tio does not change much as a function of density, since we 
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are above the critical density. The change we do see, however, 
has a temperature dependence. When densities are lowered, re- 
combination rates are lower as well. By consequence, at higher 
densities, the transition layer is closer to the edge of the cloud 
and at higher temperatures, which raises the ratio. 

In the XDR, neutral carbon occurs throughout the cloud, 
and is also abundant at relatively high temperatures. The spread 
in temperatures is large, which was already seen in Sect. 3, and 
this determines for a large part differences in the ratios. The 
temperature of a cloud is determined by Hx / n, resulting in the 
highest ratios for low n and high Fx, opposite to the situation 
in the PDR models. The three cloud types, with their different 
low, mid and high density ranges nevertheless show very simi- 
lar spreads in ratios. This is caused by the difference in column 
densities, which also has a very important effect. The low, mid, 
and high density models have their own fixed cloud size, and 
in each standard cloud type, column densities increase toward 
higher densities in the same density range. The higher density 
models contain larger regions of low temperature, which sup- 
presses the ratio at these densities even more. This can also be 
understood by considering the ratio at n = 10 4 cm~ 3 in the 
mid (type B) and high (type C) density range. The high density 
model at n = 10 4 cm -3 has a smaller cloud size and therefore 
a higher line ratio. 

5. Rotational lines 

Molecular rotational lines are also characteristic for the phys- 
ical condition of ISM gas and may also be used to constrain 
gas densities and incident radiation fields. In the following, 
we discuss a number of ratios, involving the molecular species 
12 CO, 13 CO, HCN, HNC, HCO+, SiO and CS. Although we 
reproduce in this Paper only a limited number of the diagrams 
showing the calculated line intensity ratios, all model data are 
available on-line Hence, the reader can determine all line 
ratios and integrate over all possible lines of sight as inter- 
est dictates. Here we concentrate on molecular lines that we 
consider particularly useful to distinguish between PDRs and 
XDRs. We have calculated the line intensities by using a one- 
dim ensional version of the radiation transfer code described in 
IPoelman & SpaansN2005tl2006l) . 

5. 1. CO rotational lines 

In Fig. IA.ll we show the CO(l-O) line intensity for cloud types 
A (high-density), B (mid-density) and C (low-density) for both 
the PDR and XDR model cases (cf. Table [T] All three cloud 
types are relevant for CO (and 13 CO) since these molecules are 
present ubiquitously on all galactic scales. The CO(l-O) line 
has an upper state energy E/k — 5.53 K and a critical den- 
sity n cr ~ 3 x 10 3 cm~ 3 . In the low density PDR models 
(n = 10 2 — 10 3 cm -3 ), we find that the intensity increases 
with density only. In the mid (n = 10 3 — 10 4 cm -3 ) and high 
(n = 10 4 — 10 6 5 cm~ 3 ) density range, we also find a small 
line intensity increase at higher incident radiation fields. The 
emission of the CO(l-O) line is optically thick and the emit- 

^ |http : / / www . strw . leidenuniv . nl/ ^mei jerin/grid/| 



ted line intensity mostly depends on the gas temperature near 
the r = 1 surface. Although the upper-state energy is below 
the gas temperature, we find a significant increase in the line 
intensity. At higher densities, gas-grain interactions may sig- 
nificantly increase gas temperatures in the highly attenuated 
part of the cloud. Since the escape probability j3(r) oc 1/r 
when r > 1, these parts also add a small contribution to the 
line intensity, which is proportional to ~ ln(r). For an elab- 
orate discussi on of the CO(l-O) l ine intensity dependence on 
Go and n see lWolfire et alJ dl989h . In general, the line intensi- 
ties do not vary much in the PDR models, as opposed to what 
is seen in the XDR models. In all PDR models, the CO(l-O) 
lines are optically thick. In the XDR models for type C clouds 
(n = 10 2 — 10 3 cm" 3 ), the line intensities vary over two or- 
ders of magnitude. Because of the fixed cloud size, models at 
higher densities have larger column densities. Even at the point 
farthest from the cloud edge, the low density gas in type C 
clouds causes relatively little radiation attenuation, and with 
high incident radiation it is at very high temperature and in a 
highly ionized state throughout. Under these conditions warm 
CO gas present, but only in very small amounts. Therefore, at 
such points in parameter only very weak line emission is pro- 
duced and the lines are optically thin. In the high density range 
(n = 10 4 — 10 65 cm -3 ) of type A clouds, on the other hand, 
the spread in intensity is much reduced. Although at densities 
of n = 10 4 cm~ 3 most of the cloud is still at a high temperature 
and in a highly ionized state, there is sufficient column density 
to have CO abundances large enough to produce significant line 
emission. At even higher densities (n = 10 6 cm -3 ), the col- 
umn densities are high enough to attenuate the radiation field 
in such a way that a large CO fraction is produced (~ 10 -4 ), 
but still at a temperature of T — 100 K. Here, the CO(l-O) line 
emission produced in XDRs can be two to four times stronger 
than that in PDRs. 

For the CO(2-l) line, the upper state energy is E/k — 
16.60 K and the critical density is n cr ~ 1 x 10 4 cm~ 3 . 
Although not shown, the line intensities exhibit a behavior as 
a function of density and radiation field very similar that that 
of the CO(l-O) line. Because of the higher upper state energy, 
we find a somewhat stronger dependence on radiation field in 
the PDR models. The effect of the larger critical density is hard 
to see, due to the large optical depths, but do show up when 
we consider the CO(2-1)/CO(1-0) ratio. This ratio is shown in 
Fig-S In the PDR models, the ratio does not differ more than a 
factor of two over the full density range considered here. In the 
XDR models, very large line intensity ratios of 30 or more are 
found, especially at high incident radiation fields. It is, how- 
ever, very questionable whether we will actually observe these 
high ratios, since the intensity of the emitted emission is low. 
The CO(2-1)/CO(1-0) ratio dependence on density and radia- 
tion field is in general weak, especially in PDRs, since the up- 
per state energies are not very high and the difference in critical 
density is small. 

The CO(4-3) line (see Fig. IA.2> has an upper state energy 
E/k — 55.32 K and critical density n cr ~4x 10 4 cm~ 3 . As 
expected, the emitted intensity shows more variation with den- 
sity and radiation field. In Fig.|3] we show the CO(4-3)/CO(l- 
0) line intensity ratio. In the low density range (cloud type C), 



R. Meijerink et al.: Diagnostics of irradiated dense gas in galaxy nuclei 
Low density PDRs: C0(2- 1 )/C0( 1 -0) Low density XDRs: C0(2- 1 )/CO( 1 -0) 




Density [cm 3 ] 

Mid density PDRs: C0(2-1)/C0(1 -0) 




Density [cm 3 ] 

High density PDRs: C0(2-1)/C0(1 -0) 





Density [cm 3 ] 

Mid density XDRs: C0(2- 1 )/C0(1 -0) 




Density [cm 5 ] 

High density XDRs: C0(2- 1 )/C0(1 -0) 



Density [cm 3 ] 




Density [cm ] 



Fig. 4. CO(2-1)/CO(1-0) ratio for PDR (left) and XDR (right) models. 



PDR models still produce a line ratio increase only as func- 
tion of density, but the variation in the ratio has grown to more 
than a factor of two, as opposed to only 20 percent in the cor- 
responding CO(2-1)/CO(1-0) ratio. The XDRs for this cloud 
type show a complex behavior with density and radiation field 
and the ratios cover a much larger range from about 2 to 40. At 
low radiation fields, only a density dependence is seen. At high 
radiation field strengths, the effect of the column density comes 
into play. In the mid density (type B) PDRs, the highest line ra- 
tios are seen for the highest densities and radiation fields. The 
XDRs in this range show only a dependence on radiation field. 
The effect of the higher density is compensated by the fact that 
at lower densities relatively more gas is at high temperatures. 
The gas temperature plays a large role in the density range ap- 
plicable to this cloud type. CO is present at much higher tem- 
peratures in the XDRs. Therefore, the XDR line ratios for the 
same density and incident radiation field can be more than ten 
times larger than in the PDR. This difference slowly disappears 
when the critical density of the CO(4-3) line is reached and the 



CO(4-3) line also thermalizes, which is seen in the high density 
range (cloud type A) at densities n > 10 5 cm~ 3 . 

5.2. High J -CO rotational transitions 

In cloud type A XDRs (n = 10 4 - 10 6 5 cm-' 3 ), CO is 
present throughout the cloud, even when energy deposition 
rates Hx/n are large and temperatures are high (T ~ 200 K). 
This warm CO gas produces emission originating from high 
rotational transitions even when densities are not high (e.g., 
n = 10 4 cm -3 ). Contrary to the situation in XDRs, most CO 
in PDRs is produced beyond the H/H2 transition and it has on 
average much lower temperatures (T ~ 20 — 50 K), causing 
lower intensities and line ratios. Therefore, it is very likely, that 
future missions such as Herschel/HIFI will be able to distin- 
guish between PDRs and XDRs, by observing high rotational 
transitions such as CO(16-15), CO(10-9), and CO(7-6). 

In Fig. IA.3I we show the PDR and XDR intensities of 
the CO(7-6), CO(10-9) and CO(16-15) lines for the high den- 
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Fig. 5. CO(4-3)/CO(l-0) ratio for PDR (left) and XDR (right) models 




sity range (cloud type A). We find that for both PDRs and 
XDRs, the spread in intensities increases for higher rotational 
lines, since the critical densities of these transitions are higher. 
However, this spread is much larger for PDRs than for XDRs. 
The CO(16-15) line intensity ranges from ~ 10 -10 (n = 
10 4 cm~ 3 and Go = 10 2 ) to ~ 10~ 3 erg s _1 cm~ 2 sr _1 
(n = 10 6 5 cm" 3 and G = 10 5 ) for the PDR models, while 
this is ~ 10~ 4 (n = 10 4 cm~ 3 ) to ~ 10~ 2 erg s _1 cm -2 sr _1 
(n = 10 6 - 5 cm -3 ) for the XDR models. PDRs show only sig- 
nificant CO(16-15) emission at very high densities and radia- 
tion fields (n ~ 10 6 cm -3 and Go ~ 10 4 ). This very dense and 
strongly irradiated gas, however, has a very small filling fac- 
tor on large (galaxy) scales, and, the probability of observing a 
PDR with a very high CO(16-15) intensity is low. 

The difference between the PDR and XDR models is seen 
even better by considering the intensity ratios of these high ro- 
tational transitions (see Fig.|6j. A good example is the CO(16- 
15)/CO(1-0) ratio, which ranges from 1(T 3 to 10 3 for PDRs, 
and from 10 to > 10 4 for XDRs. Especially for densities be- 



tween 10 4 — 10 5 cm 3 , it is very easy to distinguish PDRs from 
XDRs. 

5.3. 13 CO rotational lines 

The 13 CO lines have critical densities and upper state energies 
for the rotational transitions almost identical to those of 12 CO. 
As we have adopted an abundance ratio 12 C/ 13 C= 40 in our 
models, 13 CO abundances are relatively low and the lines are 
much less optically thick. 

The 13 CO(1-0) line intensities (Fig lA.4t . show the same 
trends with density and radiation field as the CO( 1 -0) line, but 
there is a larger spread in intensity. For example, the difference 
in the CO(1-0) line emission in the low density PDRs (cloud 
type C) is a factor of two and more than a factor of 3 for the 
13 CO line. This effect is even larger for the XDRs. We find sim- 
ilar results for both the 13 CO(2-l) (not shown) and 13 CO(3-2) 
line intensities as is evident from the corresponding 13 CO line 
ratios (see Fig.0and|8j. A nice illustration of the CO and 13 CO 
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Fig. 6. CO(16-15)/CO(1-0), CO(16-15)/CO(10-9), CO(10-9)/CO(7-6) and CO(7-6)/CO(3-2) ratio for PDR (left) and XDR (right) 
models. 



behavior is supplied by the the CO(4-3)/CO(l-0) versus the two and a half for the 13 CO(3-2)/ 13 CO(l-0) ratio, despite the 
13 CO(3-2)/ 13 CO(l-0) intensity ratio for the low density PDRs lower critical density of 13 CO(3-2). 
(cloud type C). The CO(4-3) line is pumped due to the fact that 
the lower rotational lines become optically thick. The CO(4- 
3)/CO(l-0) ratio changes only a factor of two, which is a factor 
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Fig. 7. i3 CO(2-l)/ i3 CO(l-0) ratio for PDR (left) and XDR (right) models. 



5.4. 13 CO/CO ratios 

In Fig. [9] and IA.6I the important isotopical intensity ratios 
of 13 CO(1-0)/CO(1-0) and 13 CO(3-2)/CO(3-2) are shown. 
Although for all cloud types (density regimes) the line inten- 
sities are larger in the XDR, the PDR isotopical ratios exceed 
those of XDRs, which means that the opacities of the XDR 
lines are larger as well. 

5.5. [CI] 609 fim/ l3 CO(2-1) ratio 

Fig.Enishows the [CI] 609 fuaJ 13 CO(2-l) ratios. For the same 
gas density and incident radiation field, PDRs have much lower 
ratios than XDRs. In the PDRs, the spread in the ratio is gen- 
erally not very large. At low densities, ratios rapidly decrease 
from 72 (n = 10 2 cm" 3 ) to 18 (n = 10 3 ) and then slowly fall 
off from 14.5 (n = 10 4 cm" 3 ) to 2 (n = 10 6 cm" 3 ). While the 
PDR ratios show a more or less steady decrease with density, a 
rather different picture is seen in XDRs. In each density range 



(cloud type), the ratio changes by several orders of magnitude. 
For the lowest density in each cloud type, the column density 
is too low to attenuate the incident radiation field sufficiently 
to allow large amounts of CO to be present. On the other hand, 
neutral carbon occurs throughout the cloud, and, therefore, a 
large increase in the ratios are seen toward high Hx / n. 

5.6. HCN rotational lines 

In Fig.ETl we show the HCN(l-O) and HCN(4-3) line inten- 
sities for high density (cloud A type) PDR and XDR mod- 
els only. Where the CO(4-3) line has a critical density of 
n cr ~ 4 x 10 4 cm" 3 , the HCN(l-O) line has a critical den- 
sity of 10 cm . Higher rotational transitions 
such as HCN(2-1), n cr ~ 4 x 10 6 cm" 3 , and HCN(4-3), 
10 cm , have even higher critical densities. The 
HCN rotational lines specifically trace the dense gas compo- 
nent in galaxies, and the line intensities from low and medium 
density gas are low. The HCN(l-O) line intensities range from 
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Fig. 8. i3 CO(3-2)/ i3 CO(l-0) ratio for PDR (left) and XDR (right) models. 



3 x 10" 11 to 2 x 10" 8 erg s" 1 cm" 2 sr" 1 for the low and mid 
density (type C en B) PDR models, and even less for the cor- 
responding XDR models. Because of their poor observational 
prospects, we have limited ourselves to only showing the model 
results in the high density range (cloud type A). Typically, the 
HCN(l-O) emission is stronger in PDRs by a factor of about 
two for densities larger than 10 5 cm" 3 . The HCN(4-3) dia- 
grams show behavior very similar to that of HCN(l-O), but the 
PDR and XDR line strengths are now somewhat closer. 

Typically, the HCN(l-O) emission is stronger in PDRs by 
a factor of about two for densities larger than 10 5 cm" 3 . Our 
results are consistent with the chemical calculations of Lepp & 
Dalgarno (1996) for different ionization rates, as follows. Our 
depth dependent models cause the HCN line emissivities to be 
the result of a line-of-sight integral over the HCN abundance 
pattern that results from a varying (attenuated) X-ray flux. Lepp 
& Dalgarno (1996, their Fig. 3) find a rather narrow range of 
ionization rates for which the HCN abundance is high and con- 
sequently the XDR HCN line emissivities have difficulty com- 



peting with the PDR ones. The HCN(4-3) contour plots show 
about the same features as seen for the HCN(l-O). However, 
the PDR and XDR line strengths are now somewhat closer. 



5. 7. HCN/CO line intensity ratios 

Intensity ratios of lines from the same species, such as the 
CO(4-3)/CO(l-0) ratio, vary with column density, due to the 
temperature gradient throughout the cloud, optical depth ef- 
fects, varying abundance etc. Intensity ratios of lines from dif- 
ferent species, in addition vary because of abundance ratio dif- 
ferences, complicating the interpretation of such line intensity 
ratios. In this section we turn our attention to the HCN/CO ra- 
tio, and start by showing in Figs. II Hand IA.8l the cumulative 
line intensity ratios for a set of PDRs and XDRs at densities 
ranging from n = 10 4 — 10 6 cm -3 and incident fluxes between 
G = 10 3 - 10 5 (F x = 1.6 - 160 erg s" 1 cm" 2 ). The cumu- 
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Fig. 9. ic, CO(l-0)/CO(l-0) ratio for PDR (left) and XDR (right) models. 



lative line intensity is the emergent intensity arising from the 
edge of the cloud to column density TVjj = n^z: 



sities Njj > 



In the XDRs, however, the variation 



between Nh = 10 and 10 cm can be rather large. 



I{z) 



1 

2^ 



K{z')dz'. 



(1) 



In the PDR, both the HCN(1-0)/CO(1-0) and the HCN(4- 
3)/CO(4-3) ratio show a minimum. The HCN abundance shows 
a drop around the H/H2 transition, while CO has its maximum 
abundance beyond the CO transition. Deeper in the cloud, the 
HCN/CO abundance ratio is more or less constant, but the CO 
line becomes optically thick and therefore the ratio increases. 

When Hx /n is low in the XDR model, the gas is molecu- 
lar, and the HCN/CO abundance ratio is more or less constant. 
A slow rise in the ratio is seen, as the CO line becomes op- 
tically thick. When the outer part of the cloud is atomic (for 
high Hx I n), HCN shows a maximum, before the H/H2 transi- 
tion. This also produces a maximum in the HCN(1-0)/CO(1-0) 
ratio. In the PDR, the variation is not that large for column den- 



In Fig. El we show the HCN(1-0)/CO(1-0) and HCN(4- 
3)/CO(4-3) line ratios for a fixed cloud size of one parsec 
(cloud type A). The variation in the ratios is relatively large, 
due to the high critical densities of the HCN transitions. The 
PDR models produce the highest HCN(1-0)/CO(1-0) ratios, 
which are attained at large densities (n > 10 6 cm~ 3 ) and may 
exceed unity. The corresponding XDRs have ratios are only 
0.1-0.2. The interpretation of the low J transitions is very diffi- 
cult due to high opacities, especially in the PDRs. For this rea- 
son, we also show the HCN(4-3)/CO(4-3) ratio, which shows 
similar trends with density and radiation field, but with some- 
what lower absolute ratios. 
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Fig. 10. CI 609 /im/ 13 CO(2-l) ratio for PDR (left) and XDR (right) models. 
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Fig. 11. Cumulative HCN(1-0)/CO(1-0) line intensity ratios for PDR (left) and XDR (right). 



5.8. HCO+ rotational lines and HCN/HCO+ line 
intensity ratios 

In Fig. \A3\ we show the HCO+(1-0) and HCO+(4-3) line 
intensities, with critical densities Tiff 2 X 10 and 4 x 
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Fig. 12. HCN(1-0)/CO(1-0) and HCN(4-3)/CO(4-3) ratios for PDR (left) and XDR (right) models. 



10 6 cm -3 , respectively. These critical densities are signifi- 
cantly lower than for HCN, causing a smaller spread in line 
intensities. Typically, the HCO + lines are stronger in XDRs 
than in PDRs by a factor of at least three. This is a direct con- 
sequence of the higher ionization degree in XDRs (Meijerink 
& Spaans 2005), leading to an enhanced HCO + formation rate. 

Note in this that Lepp & Dalgarno (1996, their Fig. 2) find 
a rather wide range of ionization rates for which the HCO + 
abundance is large. As for the HCN discussed above, this is 
consistent with our results since we integrate the depth depen- 
dent HCO + abundance profile that results from the attenuation 
of the impinging X-ray flux. The HCO + line-of-sight integral 
thus picks up a large contribution and competes favorably with 
the PDR line emissivities. 

In Figs. 1131 and IX. 101 we show the cumulative HCN(1- 
0)/HCO+(1-0) and HCN(4-3)/HCO+(4-3) line intensity ra- 
tios, for the same PDR and XDR models as in Section 5.5. 
Depending on the incident radiation field, HCN or HCO + is 
more abundant at the PDR edge of the cloud. Around the 
H/H2 transition a minimum in the HCO + abundance is seen 
in the PDR. Deeper into the cloud, the HCO + abundance 
increases again, and is then constant. This is also the case 
for HCN, and the HCN/HCO+ abundance ratio is larger than 
unity. Therefore, at sufficiently large columns and densities, the 
HCN(1-0)/HCO + (1-0) line intensity ratio becomes larger than 
one. 

In the XDR models, HCO + is chemically less abundant 
than HCN for very large Hx/n (Meijerink & Spaans 2005). 
For larger columns HCO + becomes more abundant, however, 
and eventually the cumulative column density of HCO + be- 
comes larger than HCN (see specifically Fig. 10 in Paper I). 



This follows directly from the fact that the HCO + abundance 
is high over a much wider range of ionization rates than HCN 
(Lepp & Dalgarno 1996, their Figs. 2 and 3). 

Fig- El clearly shows that the HCN/HCO+ ratio discrim- 
inates between PDRs and XDRs in the density range be- 
tween n = 10 5 and 10 6 5 cm" 3 (cloud type A). The HCN(1- 
0)/HCO+(1-0) and HCN(4-3)/HCO+(4-3) line ratios are both 
much larger in the PDR models than XDR models, for columns 
of 10 23 cm~ 2 and larger (Paper I). The difference ranges be- 
tween a factor of 4-10, depending on the density. The XDR 
HCN(1-0)/HCO+(1-0) ratio becomes larger than unity for 
more modest columns of 10 22 5 cm~ 2 and less. 



5.9. HNC/HCN ratios 

The critical densities of HCN and HNC are almost identical, 
so that the only differences in line ratio should be due to dif- 
ferences in the abundances. In PDRs, HCN is more abundant 
in the radical region, but deeper in the cloud the abundance ra- 
tio approaches unity. In XDRs, HCN is more abundant in the 
highly ionized part of the cloud. However, HNC is equally or 
even more abundant than HCN deep into the cloud. As a result, 
the HNC(l-0)/HCN(l-0) line intensity ratio is around one for 
the PDRs if the column density is larger than 10 22 cm -2 , while 
the HNC( 1 -0)/HCN(l-0) ratio is less than unity for N H < 10 22 
cm -2 . The XDR models, however, show low ratios for the low, 
~ 10 4 cm~ 3 , densities and strong, > 10 erg s _1 cm~ 2 , ra- 
diation fields. The ratios increase for lower incident radiation 



fields, and at highest densities (n = 10 cm 3 ) the line ratios 
are always larger than one, irrespective of irradiation. 
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Fig. 13. Cumulative HCN(1-0)/HCO+(1-0) line intensity ratios for PDR (left) and XDR(right). 
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Fig. 14. HCN(1-0)/HCO+(1-0) and HCN(4-3)/HCO+(4-3) ratios for PDR (left) and XDR (right) models. 
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Fig. 15. Cumulative HNC(l-0)/HCN(l-0) line intensity ratios for PDR (left) and XDR(right). 
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Fig. 16. HNC(l-0)/HCN and HNC(4-3)/HCN(4-3) ratios for PDR (left) and XDR (right) models. 



In the PDRs, the HNC(4-3)/HCN(4-3) ratio quickly drops 
below unity at densities below n = 10 5 cm" 3 . This density 
is far below the critical densities of the lines, and therefore 
high temperatures are needed to excite them. Such high tem- 
peratures are indeed found in the radical regions of the PDRs, 
but there the HNC abundance is much lower than the HCN 
abundance, which explains the drop in the ratio. In the XDRs, 
the HNC(4-3)/HCN(4-3) ratios are quite similar to the HNC(1- 
0)/HCN(l-0) ratios, except for densities n > 10 6 cm -3 , where 
they are even high than these. 

5.10. SiOandCS 

Although SiO is usually considered to be a good tracer of 
shocks, we do find that the SiO(l-0)/CO(l-0) ratio is typically 
larger in XDRs than in PDRs by a factor of 2-3, for densi- 
ties around 10 5 5 cm" 3 . For the higher excitation J=4-3 lines, 
the effect disappears because CO is generally warmer in XDRs 
compared to PDRs. 

The CS(l-0)/HCN(l-0) ratio is a factor of two larger 
(smaller) in XDRs for densities above (below) 10 5 cm" 3 . 
Interestingly, the corresponding 4-3 ratio continues this trend 
but changes in the ratio from 10 4 to 10 6 cm" 3 are now as large 
as a factor of 10. 

6. Column density ratios 

Unfortunately, for many molecular species of interest no reli- 
able collisional cross sections are available. For these species 
we are unable to accurately predict line intensities, but we can 
still calculate the column density ratios. In this section, we dis- 



cuss column density ratios for a number of species that are of 
potential interest in attempts to discriminate between PDRs and 
XDRs. 

6. 1. CN/HCN column density ratio 

In Fig.^J we show the CN/HCN column density ratios. There 
is an enormous difference between the ratios for PDRs and 
XDRs. In the PDRs, the cloud type A ratios range from 0.5 
(n ~ 10 6 cm" 3 ) to 2.0 (n - 10 4 cm" 3 ), while in the XDR 
models the same ratio varies from 40 (n ~ 10 6 cm -3 ) to over 
a 1000 (n ~ 10 4 cm" 3 ). We find higher CN/HCN ratios at 
lower densities, where the chemical rates are lower, making it 
more difficult to form large molecules. The PDR ratios are only 
dependent on density, which is explained by the fact that most 
CN and HCN molecules are formed beyond the H/H2 transi- 
tion. This part of the cloud is shielded from FUV photons and 
the chemistry is dominated by the cosmic ray ionization rate, 
which is the same in every model. In the XDR models, there is 
much less variation in the CN and HCN abundance throughout 
the cloud. The variations do not exceed more than two to three 
orders of magnitude, while this is over ten orders of magnitude 
in the PDR models. All parts of the cloud contribute almost 
equally to the column density ratio, including the region with 
very high Hx /n>. Hx /n is a major factor in the resulting ra- 
tio, and therefore the XDRs also show a large dependence on 
incident radiation field. 

On the website, we show the PDR and XDR cumulative 
CN/HCN ratios for a few specific densities and radiation fields. 
The variation in the cumulative column density ratio is much 
less in the XDR models than in the PDR models. At PDR 
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Fig. 17. CS(l-0)/HCN and CS(4-3)/HCN(4-3) ratios for PDR (left) and XDR (right) models. 
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Fig. 18. CN/HCN column density ratios for PDR (left) and XDR (right) models. 



edges, the gas is highly ionized (as in the XDRs), and here we 
find ratios resembling those of XDRs. Abundances, however, 
are very low here because of the high photo-dissociation rate. 



6.2. CH/HCN column density ratios 



In Fig. IA.15l we show the CH/HCN column density ratio. The 
differences between the ratios in PDRs and XDRs are even 
larger than in the case of CN/HCN. While the PDR ratios in- 
crease from 0.2 (n ~ 10 6 cm~ 3 ) to 0.9 (n ~ 10 4 cm -3 ), the 
XDR ratios range from 20 to more than 10000. The PDR ra- 
tio does not depend on density only. At relatively low densi- 
ties (n = 1 — 3 x 10 4 cm -3 ), we find a dependency on in- 
cident radiation field as well. CH reaches its maximum abun- 
dance at greater depths than HCN, which means that the cu- 
mulative ratio is still increasing at large column densities. Our 



model cloud sizes are not large enough to allow the ratio to 
converge to a constant ratio (see website). For the XDR mod- 
els, we can roughly state that the ratio increases toward higher 
Hx /n. However, at n ~ 10 5 cm -3 the lowest ratio is not found 
at the lowest incident radiation field strength. At such densities, 
the limited cloud sizes are comparable to the depth of the H/H2 
transition. This transition is included in models with the low- 
est radiation field strengths, but not in those with the highest 
radiation field strengths. 



6.3. Chi^/HCN column density ratios 

In Fig. lAT6l we show the CH+/HCN ratios for both PDR 
and XDR models. They range from 10 -6 (n ~ 10 6 cm -3 ) 
to 6 x 10 -4 (n ~ 10 



4 cm" 3 ) in the PDRs, and from 10 -3 

-3\ 



(n ~ 10 b cm" 3 ) to over 1000 (n ~ 10 4 cm" 3 ) in the XDRs 



R. Meijerink et al.: Diagnostics of irradiated dense gas in galaxy nuclei 



19 



In the PDR models, the highest CH + abundance is seen close 
to the edge of the cloud but it decreases very quickly beyond 
the H/H 2 transition. HCN, on the contrary, reaches its highest 
abundance beyond this transition. This explains the decrease 
of the PDR cumulative column density ratios with increasing 
depth (see website). The PDR models show a dependence on 
both density and incident radiation field for this ratio. The de- 
crease in the ratio with increasing density is caused by the 
fact that on the one hand CH + is more easily destroyed (due 
to higher recombination rates) and on the other hand HCN 
more easily formed when densities are higher. The ratios do 
not converge to a constant value deep into the cloud. At identi- 
cal densities, the ratio increases with increasing radiation field, 
since the HCN column densities become smaller while the 
CH + column densities become larger for the cloud-size con- 
sidered here. In the XDR models, the largest ratios are seen 
for the highest Hx/n (low density and high incident radia- 
tion field). The fluctuations in the CH + and HCN abundances 
are more gradual. The HCN abundance increases and the CH + 
abundance decreases when the X-ray photons are gradually ab- 
sorbed. Therefore, the XDR cumulative column densities ratios 
show less variation than the PDR ratios (see website). 

6.4. HCO/HCO+ column density ratios 

Fig. IA. 171 shows the HCO/HCO + column density ratios. For 
the PDRs, we find much larger ratios than for the XDRs. The 
PDRs show ratios between 0.1 (n ~ 10 4 cm -3 ) and 10.0 
(n ~ 10 6 cm -3 ), while the ratios in the XDR range from 10 -5 
(n ~ 10 4 cm -3 ) to 10 -3 (n ~ 10 6 cm" 3 ). In the PDRs, the 
ratios depend only on density for radiation fields Go < 10 4 . 
With larger incident radiation fields, the ratio becomes depen- 
dent on the radiation field strength as well. The HCO abun- 
dance reaches its maximum and more or less constant abun- 
dance somewhat deeper into the cloud than HCO + . As col- 
umn density increases, the HCO/HCO+ ratio slowly converges 
to a constant value (see website). For large radiation fields 
(Go > 10 4 ), the cloud size considered here is too small to 
allow column densities to converge to a constant ratio. In the 
XDR models, the ratio depends on both density and radiation 
field in all regimes. The lowest ratios are seen for high Hx/n 
(n = 10 4 cm -3 and F x = 160). 

6.5. HOC+/HCO+ column density ratios 

In Fig. IA.18I we show the HOC + /HCO + column density ra- 
tios. The XDRs show larger ratios, ranging from 10 -4 (n ~ 
10 6 cm" 3 ) to 0.6 (n ~ 10 4 cm -3 ), than PDRs, where the ratios 
range from 10 -7 to 10 -5 . The XDR ratios increase for larger 
Hx/n (with maximum at n ~ 10 4 cm -3 and Fx = 160). 
The behavior is more complex for the PDR models. The ratios 
depend on density only for Go > 10 4 , but below this value 
there is also a dependency on radiation field strength. This is 
explained by the fact that ratios drop very fast, when the gas 
becomes molecular (see website), but abundances are still sig- 
nificant. A large fraction of the gas at the edge of the cloud is 
molecular for densities n > 10 5 5 cm -3 and radiation fields 



Go < 10 4 . Here we find a fast decrease in the ratio for lower 
radiation fields. 

6.6. NO/CO column density ratios 

In Fig lA.19l we show the NO/CO column density ratios. The 
ratios in XDRs are much larger, 10 -4 (n ~ 10 6 cm -3 ) to 10 -3 
(n ~ 10 4 cm- 3 ), than for PDRs (10 -6 - 10" 5 ). In the PDR 
the ratios are largely determined in the molecular part of the 
cloud and depend more or less on density only. The ratios in 
the XDR depend mostly on Hx /n. The largest ratios are seen 
for the largest H x /n (n ~ 10 4 cm" 3 and Fx = 160). 

CO(1-0) lines have optical depths of r(CO(l - 0)) - 100. 
Therefore, it could be possible to observe NO(1-0)/CO(1-0) 
line intensity ratios as large as 0. 1 , while the maximum column 
density ratios are only 10 -3 . 

6.7. N 2 hl + /CO column density ratios 

In Fig. IA.20l we show the N 2 H + /CO column density ratio. In 
both PDRs and XDRs the model ratios are rather low. They de- 
crease with density in the PDR (10 -7 — 10 -9 ), which is oppo- 
site to the behavior in XDRs. In the XDRs, we find a decrease 
with increasing H x / n. There are collisional data available for 
N 2 H + , but since the abundances are so low, the line intensities 
are too small to be observable. 

7. Summary and outlook 

We have presented a large set of PDR and XDR models that 
can be used to determine the physical conditions that pertain to 
irradiated gas clouds. This grid spans a large range in densities 
(ra H = 10 2 - 10 6 - 5 cm -3 ), irradiation (G = 10 5 - 10 5 and 
Fx = 1.6 x 10 -2 — 160 erg cm -2 s -1 ) and column densities 
(A H = 1.5 x 10 22 - 1 x 10 25 cm -2 ). We have used the results 
to make predictions for the intensities and ratios of the most 
important atomic fine-structure lines, e.g., [CII], [OI], [CI], 
[Sill], and [Fell], rotational lines for molecular species such as 
HCO+, HCN, HNC, CS and SiO (up to J = 4), CO and 13 CO 
up to J = 16, and for column densities for CN, CH, CH + , 
HCO, HOC+ , NO, and N 2 H+ . It is not possible to to present all 
the results, but they are available on-line at the following URL: 
|http : / /www . strw . leidenuniv . nl/~mei jerin/grid/| 
Here we summarize the most important conclusions: 

1 . The surface temperatures are higher (lower) in PDRs com- 
pared to XDRs for densities n > 10 4 cm -3 (n < 
10 cm -3 ). Two opposing effects play a major role in de- 
termining the resulting surface temperature: (1) The heat- 
ing efficiency, which is much higher in XDRs (up to 70 
percent) than in PDRs (0.5-3.0 percent); (2) The absorption 
cross sections which are much smaller for X-rays than for 
FUV photons. 

2. For the atomic lines, we find that the fine-structure line 
ratios of [Sill] 35 /jm/[CR] 158 /im, [OI] 63 /im/[CII] 
158 ^m, [Fell] 26 um/[Cll] 158 /im, and [CI] 369 um/[Cl] 
609 um are higher in XDRs than in PDRs, for a given den- 
sity, column, and irradiation strength. Whereas PDR ratios 
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depend on density and irradiation strength only, XDR de- 
pend on column density as well. In PDRs, fine-structure 
line emission is only produced at the edge of the cloud, 
while in XDRs almost all parts of the cloud contribute. 

3. We find higher CO line ratios for XDRs. In PDRs, CO is 
formed beyond the H/H2 transition and typically has tem- 
peratures in the range T ~ 20 - 50 K. In XDRs, CO 
is present throughout the cloud in significant abundances, 
even in the highly ionized part. When using CO line ratios, 
the best way to distinguish between PDRs and XDRs is to 
consider ratios such as CO(16-15)/CO(1-0), where the dif- 
ferences are largest (Fig. [6}. 

4. HCN/HCO" 1 " ratios discriminate well between PDRs and 
XDRs, even in the lower rotational lines, and especially 
when densities are as high as n > 10 5 cm -3 . At such densi- 
ties, the HCN(1-0)/HCO + (1-0) line intensity ratios are < 1 
in XDRs, while PDRs have ratios > 1 for column densities 
iV H > 10 23 cm" 2 . Although the HCN/HCO+ line ratio 
in an XDR may become even larger in clouds of relatively 
modest density (10 4 cm -3 ) subjected to high radiation field 
strengths (> 100 erg s -1 cm -2 ), we find that the line in- 
tensities in this part of parameter space are too low to be 
detectable (see also Meijerink et al. 20061). 

5. For densities between 10 4 and 10 6 5 cm -3 , HNC(1- 
0)/HCN(l-0) ratios in PDRs are of order 1 « 1) for 
columns larger (smaller) than 10 22 cm -2 , while the ratios 
range between 0.2-1.2 for XDRs. PDR HNC(4-3)/HCN(4- 
3) ratios never exceed unity, while we find XDR ratios up 
to 1.6 for densities n > 10 6 cm -3 . 

6. HCN/CO ratios are typically smaller for in XDRs than in 
PDRs, for two reasons: (1) The HCN abundance is boosted 
only in high (column) density gas, with columns in excess 
of 10 23 cm -2 and densities larger than 10 4 cm -2 ; (2) CO 
is warmer in XDRs, which leads to stronger emission, and 
this suppresses the HCN/CO ratios as well. In PDRs, the 
very high HCN(1-0)/CO(1-0) ratios of order > 1 are only 
obtained at very high column densities (JVh > 10 23 cm -2 ). 
High column densities Njj are needed to obtain a large 
HCN column density and to obtain optically thick HCN 
lines. This happens much faster for CO, due to its higher 
fractional abundance. 

7. We find that CN/HCN, NO/CO, and HOC+/HCO+ column 
density ratios are discriminant between PDRs and XDRs. 
Molecules such as CH and CH + also look very promis- 
ing. However, we urgently need reliable collisional cross 
section data in order to make predictions for observed line 
ratios. 

We conclude that both atomic fine-structure and molecular ro- 
tational lines have significant diagnostic value to allow us to 
distinguish between clouds irradiated by star-bursts (FUV) and 
by active galactic nuclei (X-ray) in the centers of galaxies. 
Note, however, that the line intensities (and ratios) presented 
in this work do not take the, likely complex, kinematics of nu- 
clear gas into account. For example, gas rotating in an accretion 
disk around an AGN will cause the line widths of, say, HCN, 
HNC and HCO + transitions to differ, depending on where their 
chemical abundances peak with depth. 



The XDR/ AGN contribution will typically be of a much 
smaller (possibly beam diluted) angular scale than that 
of a PDR/Starburst. A 10-25% PDR contribution may al- 
ready suppress our ability to recognize XDR excitation from 
HCN/HCO+ and HNC/HCN line ratios. A solution to this can 
be found in the very high J CO lines (e.g., CO(16-15)), that 
are excellent indicators of an XDR contribution. These very 
high rotational lines will be observable with the ESA Herschel 
(HIFI) space observatory scheduled for launch in the near fu- 
ture; they can be seen in absorption in the near-infrared with 
Subaru. Currently available (sub)millimeter facilities lack the 
spatial resolution to separate the PDR/stellar and XDR/ AGN 
contributions in distant active galaxies. Up to now, these com- 
ponents can only be spatially resolved in the Milky Way. 
However, the resolving power of ALMA will bring this pos- 
sibility within reach for external galaxy centers as well. 
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Fig. A.l. CO(1-0) intensity for PDR (left) and XDR (right) models. 
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Fig. A.2. CO(4-3) intensity for PDR (left) and XDR (right) models. 
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Fig. A.3. CO(16-15), CO(10-9), and CO(7-6) intensity for PDR (left) and XDR (right) models. 
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Fig. A.4. 13 CO(1-0) intensity for PDR (left) and XDR (right) models. 



26 



R. Meijerink et al.: Diagnostics of irradiated dense gas in galaxy nuclei 

Low density PDRs: ,3 C0(3-2) intensity Low density XDRs: ,3 C0(3-2) intensity 

1.0 : 




100 



Density [cm 3 ] 

Mid density PDRs: 13 C0(5-2) intensity 



1000 



1000 



















s 

\ ; 
V 



Density [cm 3 ] 

High density PDRs: ,3 CO(3-2) intensity 





100 



Density [cm 3 ] 

Mid density XDRs: ' 3 C0(5-2) intensity 



1000 




E 

o 

10 



Density [cm 3 ] 

High density XDRs: l3 C0(5-2) intensity 



] o 

/ 


a? '* ^ 


b \ 
x \ 


K : 


o \ 




■ / \ . \ 





10 5 

Density [cm -3 ] 



10" 



10 s 

Density [cm -3 ] 



10 6 



Fig. A.5. 13 CO(3-2) intensity for PDR (left) and XDR (right) models. 
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Fig. A.6. 13 CO(3-2)/CO(3-2) ratio for PDR (left) and XDR (right) models. 
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Fig. A.7. HCN(l-O) and HCN(4-3) intensity for PDR (left) and XDR (right) models. 
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Fig. A.8. Cumulative HCN(4-3)/CO(4-3) line intensity ratios for PDR (left) and XDR(right). 
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Fig. A.9. HCO+(1-0) and HCO+(4-3) intensity for PDR (left) and XDR (right) models. 
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Fig. A.10. Cumulative HCN(4-3)/HCO+(4-3) line intensity ratios for PDR (left) and XDR(right). 
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Fig. A.ll. Cumulative HNC(4-3)/HCN(4-3) line intensity ratios for PDR (left) and XDR(right). 
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Fig. A.12. SiO(l-O) and SiO(4-3) intensity for PDR (left) and XDR (right) models. 
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Fig. A.13. SiO(l-0)/CO and SiO(4-3)/CO(4-3) ratios for PDR (left) and XDR (right) models. 
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Fig. A.14. SiO(l-O) and SiO(4-3) intensity for PDR (left) and XDR (right) models. 
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Fig. A.15. CH/HCN column density ratios for PDR (left) and XDR (right) models. 
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Fig. A.16. CH+/HCN column density ratios for PDR (left) and XDR (right) models. 

PDRs: HC0/HC0 + column density ratio XDRs: HC0/HC0 + column density ratio 




Fig. A.17. HOC/HCO+ column density ratios for PDR (left) and XDR (right) models. 

PDRs: H0C + /HC0 + column density ratio XDRs: H0C + /HC0 + column density ratio 




Density [cm 3 ] Density [cm 3 ] 

Fig. A.18. HOC+/HCO+ column density ratios for PDR (left) and XDR (right) models. 

PDRs: NO/CO column density ratio XDRs: NO/CO column density ratio 
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Fig. A.19. NO/CO column density ratios for PDR (left) and XDR (right) models. 




Fig. A.20. N 2 H+/CO column density ratios for PDR (left) and XDR (right) models. 



